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OBJECTIVE: 


An  unprecedented  amount  of  Lagrangian  data  is  now  being  collected  by 
oceanographers.  The  general  feeling  is  that  such  data  are  important  for  an  accurate 
description  of  ocean  circulation.  Regrettably,  little  hard  analysis  has  gone  into 
deployment  strategies  and  even  less  into  quantitative  use  of  Lagrangian  data  after  it 
has  been  collected.  This  is  unfortunate,  since  such  data  can  enhance  numerical 
modeling  efforts  through  data  assimilation,  which  makes  explicit  use  of  the  data’s 
Lagrangian  character.  These  data  Cem  also  be  used  to  assess  a  model’s  ability  to 
predict  the  evolution  of  transient  events. 

The  long  range  objective  of  this  project  is  to  develop  theoretically  based 
quantitative  procedures  for  utilizing  Lagrangian  data  for  both  of  these  applications. 
The  focus  of  this  phase  of  the  effort  will  be  on  using  both  surface  and  subsurface 
Lagrangian  data  for  a  quantitative  assessment  of  circulation  models. 

APPROACH: 


The  approach  used  here  is  an  improvement  of  that  used  by  Eremeev  et  al.  {Jour. 
Geophys.  Res.  97:  9733-9753,  1992).  This  is  summarized  briefly  here. 

The  basis  for  the  procedure  is  a  general  Eulerian  representation  of  any 
hydrodynamic  field  P’’  as: 

r'(x,,,i)  =  f;>t;(,,f)r;(x,t)  (i) 

il;=l 


Here,  the  superscript  r  could  represent  components  of  the  velocity  field  and  77  is  a 
surface  that  is  a  chciracteristic  of  either  a  model  or  the  Lagrangian  observations. 
This  might  be  depth  for  a  numerical  model,  and  an  isopycnal  or  isobaric  surface  for 
data.  The  Lagrangian  representation  of  the  same  field  simply  replaces  the  horizontal 
position  X  with  the  path  data  for  Lagrangian  probe  p. 

The  key  issue  is  the  specification  of  the  basis  function  P]^  in  (1).  A  number  of 
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options  axe  available,  but  the  one  most  appropriate  for  this  effort  is  to  utilize  a 
representation  of  the  velocity  field  as: 

u  =  V  X  [nV*  +  V  X  (n<^)l  (2) 

where  n  is  normal  to  t].  Projection  of  u  and  V  ^  u  onto  rj  give: 


n  •  u  =  -  V  [(n  •  n)  V  —  n  (n  •  v^)] 

n  •  (v  X  u)  =  n  •  Vh  (n  X  v^)  -  V  ((»  •  n)  V  V'  -  n  (n  •  vV’)]  (3) 

The  solution  of  (3)  involves  eigenfunctions  of  a  Poisson  problem  which  ensure  no 
slip  and  no  normal  flow  at  the  domain  boundaries  -  typical  numerical  modeling 
constraints.  These  eigenfunctions  are  the  basis  functions  Tl  in  (1).  They  can  be 
calculated  (with  some  difficulty  and  care)  for  any  given  model  geometry.  Once 
determined,  they  are  used  for  both  model  output  and  observations. 

PROGRESS  SUMMARY: 

We  have  been  working  on  this  project  for  less  than  6  months,  so  only  preliminary 
results  axe  available.  Although  the  approach  described  above  is  well  documented  in 
the  literature,  there  is  precious  little  experience  with  applications.  In  order  to 
develop  some  focus,  we  elected  to  first  apply  the  method  to  enclosed  or  nearly 
enclosed  basins  with  a  substantial  amount  of  Lagrangian  data.  We  have  succeeded 
so  far  in  the  first  task  of  developing  numerical  methods  for  calculating  the 
eigenfunctions  in  (1)  above  for  no  slip  boundary  conditions. 

APPLICATIONS: 

When  developed,  this  technology  will  impact  naval  operations  in  two  ways.  First,  it 
will  provide  a  powerful  quantitative  tool  for  comparing  circulation  model  predictions 
with  Lagrangian  observations.  This  would  assist  in  the  interpretation  of  data  from 
transient  events  and  provide  a  more  powerful  method  for  assessing  model 
performance  than  is  currently  used.  Second,  it  will  provide  a  means  for  assimilating 
drifter  and  float  data  into  predictive  models.  Such  a  development  would  correct  a 
long-standing  problem  with  present  assimilation  methods,  which  don’t  take  full 
advantage  of  the  Lagrangian  character  of  the  data. 
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